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Abstract 


In this article, the chemical kinetics problem is presented. The Differential Transformation Method (DTM) is presented 
to find an approximate analytical solution to the chemical kinetics problems. Proposed problem is system of nonlinear 
ordinary differential equations. We show converges analysis of the differential transform method. DTM applied to two 
examples. The method shows the form of fast converging series and the results prove the applicability of the proposed 
method, which gives accurate results. 
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1 | Introduction 


In this paper we are concerned chemical kinetics problems system; 
COD iconsce 
Computational dx, 
Algorithms and dt. i a a 
Numerical Dimensions. dx» 2 1 
This article is an open “dt. le a a a @) 
access article distributed dx 3 2 
under the terms and dt = gee 
conditions of the Creative 
Commons Attribution With the initial conditions; 
(CC BY) license 
(http://creativecommons. at 0) =1, m2 0) =O, x3 0) = 0. (2) 
org/licenses/by/4.0). Whete a4, a> and ay are the reaction rates. In chemistry the chemical kinetics system is well described 


by nonlinear system of ordinary differential equations. Robertson [4] introduced the system in 1966, 
indicating the mathematical model of a chemical kinetic problem [1], [2]. 
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In recent years, numerous works have been focusing on the development of more advanced and efficient 
methods for the Homotopy Perturbation Method (HPM) and the Variation Iteration Method (VIM) [2], 
the picard-pade technique [3], the homotopy analysis method [3] and the Adomian’s decomposition 
method [5] and iterative method [6]. The Differential Transformation Method (DTM) was first created by 
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Zhou [7] to obtain approximate-analytical solutions to ordinary differential equations using Taylor series 
formulation. The differential transform method obtains an analytical solution in the form of a polynomial. 
It is different from the traditional high order Taylor’s series method, which requires symbolic competition 
of the necessary derivatives of the data functions. The results of applying DTM to the chemical kinetics 
problems will be presented. 


2 | Differential Transform Method 


Differential transform of function y(x) is defined as follows: 


dky(x) 
dxk 


1 
YK=5 (3) 


In Fg. (3), y x)is the original function and Y(k) is the transformed function, which is called the T-function. 
Differential inverse transform of Y(k) is defined as 


co 


yx)= y xkY k). (4) 


k=0 
Prom Eq. (3) and (4), we obtain 


dky(x) 
yx = ae dxk | ©) 


Eq. (2) implies that the concept of differential transform is derived from Taylor series expansion, but the 


method does not evaluate the derivatives symbolically. 


However, relative derivatives are calculated by an iterative way which are described by the transformed 
equations of the original functions. In this study we use the lower case letter to represent the original 
function and upper case letter represent the transformed function. 


From the definitions of Egs. (2) and (4), it is easily proven that the transformed functions comply with the 
basic mathematics operations shown in Tabk 7. 


In actual applications, the function y(x) is expressed by a finite series and Eq. (4) can be written as: 
m 


yx)= > xy k). (6) 
k=0 
Eq. (5) implies that Yeni“ Y 4), is negligibly small. In fact, m is decided by the convergence of natural 
frequency in this study. 


Table 1. The fundamental operations of DTM. 


Original Function Transformed Function 


y(x) = u(x) + v(x) Y(k) = U(k) + V(K) 
y(x) = cw(x) Y(k) = cW(k) 
y(x) = dw/dx Y(k) = (K + 1)W(k + 1) 


y(x) = diw/dxi = -Y(k) = (K +1)(K + 2)... K+j)W(k + j) 
y(x) = u(x)v(x) is 
Y(k) = >> U(r) V(K - 1) 
r=0 


Theorem 1. If y x) = exp(x) then k) =— 


Proof. By using Eg. (3), we have 


1 d(e* 1 1 
Wk)=- al =—e*] =—, 
Theorem 2. If y x) = 2 ~~ then Y k) =“ k +m). 
Proof. By using Eg. (3), we have; 
1 dk o™ux 1o™**a x k+m 
Yk)= k = k ah Uk+m) 
kKidx® dx™ Jug ki Oxe* | 4 k! 
— ym = _ = 1 k=m 
Theorem 3. If y x) =x" then k) =6k m)= {i cen 
Proof. By using Eg. (3), we have; 
1 0K(xK) oy 
We <2?) oleae tet) k=m 
kK! axkK hog) 1 ak x™) _ k#ém 
kl axk 7” 
wk 
Theorem 4. If w x) = sin(wx + a) then W k) = ~ sin(= +a). 
Proof. By using Eg. (3), we have; 
1 0 sin wx+a) 1 1 : Tl 
k=1 W1)= qo ae |p = 7, weos wx + a = Ewsin(( + cr) + wx} - 
1 : [= ‘i ), 
=—wsin|—+a 
1! 2 ; 
1 0*sin wx+a 1 
k=2 W2)= i = 5,w? cos (wx + +0] 
Ox -0 2! x=0 
1 2 a 1 2... (2% 1 2 
= sin(5 + fe + cr) + wx)| = sin(> +at wx)| = 5. sin (= + a). 


In the general form we have; 


1a” + 
k=n Wn= —e ») 


1 n.: nim 1 
=—Ww sin(wx + = + a)| =—w *sin(S +a). 
x=0 =™ x=0 


Theorem 5. If w x) = cos(wx + a) then W k) = © cost +a). 


Proof. By using Eg. (3), we have; 


noses wx + a) -1 
———|  =—wsin wx+a) 


k=1 W1 
i Ox say | ee, 
iw con(5 + “ naw +8) 
—w cos a}|+wx =—wWw COS al, 
~ 7 vey 2 
1 0%cos wx +a = Tl 
k=2 W2)= : _ w? sin (wx +5 +a] 
2! Ox x0 2! x=0 
ae T at Bane (S a 1 Bane (2 
=5w cos(+ + (5 + or) +wx)| =a cos (5 + a + wx) Lae rl cos (= + ox) 
In the general form we have; 
1 woes wx+a) log nm == 
k=n Wno)= ie ae =a cos(wx + 7 +o)| /w "cos(F +c). 


In this part, we will show converges of the differential transform method to the exact solution by proof 


of some theorems [8]. 


Theorem 6. Consider the first order ordinary differential equation in here m,n are constants; 


dw(x) 
m at nw x) = 0,w 0)=y. (7) 
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The differential transform method converges to the exact solution. 
Proof. The general solution of differential equation is as w x) = Be™/", 


Now by taking the differential transform of Eg. (7), we have; 


m W[K + 1] +nW[K] = 0, W[0] = 6. 


—nW[K] ‘ : 1 k 
We have W[K +1] = =f By solve this recurrence relation, we get W[K] = i (=) B. 


By using Eq. (8), the solution obtained by DTM is as follows: 


1 /-nx\k ees 
w x) = 8 Dog (=) =fe m, (8) 
Theorem 7. Consider the homogeneous ordinary differential equation of order n, where v;,i =1,2,...,n 
ate variables; 


Vo X)w+ Yea Vr x) = 0. (?) 
With initial conditions 


diw 
Ww x)= Bo, “aa 


=,i=1,2,..,n-1. (10) 
x=0 


The differential transform method converges to the exact solution. 


Proof. By using with both Tabé 7 and the Eq. (9), we have; 
Dk -o (VolkIWEk =m] + D2, [VelKI[Ifa km +p) ]Wik-m +r]]) =o. 
From the initial conditions, we have; 
WI0] = fo, WIi] = 2,1 =12,....n-1. 
By substituting the differential transforms, we have; 
w x) = Deo WIk]xk, v, x) = No A rl k]x*. 
Substituting the above formulas in Eg. (9), we get; 
(oe) oO [oe] d* [oe] 
0 = Dido A olkIx* Dido WIkIx* + DE (xz A Ak}x*. (5 Deo WIkIx)) 
= Do Daneo A olm] Wik = m)x* D2 (Deo A clk]xk. D9 Wim + r])(TIEo m + 
r—}))x™) = Do DK-o(A olm]WIk - m] + D2, A -[m]W[k - m +r], k- 
m+ i) )xk. 
We comparing coefficients and get; 


yk _) A ofm]W[k — m] + 52, A ,[m]W[k - m+ rT k-m+j)=0. 


3 | Numerical Example 


Example 1. We consider the system of chemical kinetics; 


dx, 7 
—— = —a1X1 + a9X0X3, 
dt 1X1 2X2X3 
dx, _ , (11) 
Ty, — 43X14 — 44X2X3—A5XQ", 
iat 
X 
“a 7 BD 
With the initial conditions 
x, 0) =1,xz 0) =0, x3 0) =0. (12) 


When taking the one dimensional differential transform of Eg. (17), we can obtain: 


X,;n+1)= went a,X1 n) +a. Yio Xo r)X3 n-1)], 
Xgnt+laq anl@% N) — a4 Yr-o Xz 1)X3 N- 41) — as Yr-oX2 1X2 n—1)], (13) 
X3n+1)= ma [ag Uir-0 X2 r)X2 n-1)], 
The initial conditions can be transformed at xg = 0, as; 
Ay OS 1' yp Xo OH 0, Xe) HO. (14) 
Porn = 0,1,2,..., X,(1), Xp(n), X3(1), coefficients can be calculated from Egs. (13) and (14); 


X, 1) =-ay, Xp 1) =a3, X31) =0, Xy 2) => 


(15) 
X2 2) =—S, X3 2) =0, X13) == 


Xo a= (ae Lae —as a3) *), X3 3) = —— as) yee & 


By using the inverse transformation rule for one dimensional in Eg. (6), the following solution can be 


obtained: 
aio aris 
X, t) =l-ajt+— t*-—t°t 
1 ) I 2 PS 
a3a4 1 (ag ay)? 
Xo t) =a3t- 5 eae 5 —as5 as)? tte, (16) 


2 
ie po 


Example 2. We consider the system of chemical kinetics; 


dx, 
at ' 
se 2 17 
Sse, (17) 
dt 1 A2 
dx3 = 2 
dt 
Subject to the initial conditions; 
xX, 0) =1,x2 0) = 0,x3 0) = 0. (18) 
When taking the one dimensional differential transform of Eqs. (17) and (18), we can obtain: 
X;n+1)= wen X, n)], 
n — 
X, n+) =z i aap PX 0) - Dy Xe DX. n=), as 
X3n+1)= =) Xp r)X2 n- of 
And 


Xy DHT Xs 0) =0,X3 0) =0. (20) 
Porn = 0,1,2,..., X,(n), X2(n), X3(n), coefficients can be calculated from Egs. (19) and (20). 
x, 1) = -1,x2 1) = 1,x3 
1 =l 
Xi 2)=5 7X2 2) => 1X3 2) = (21) 
= eG 1 
X1 3) =—,X2 3) = Bx 3) = 


By using - inverse ae acd a for one dimensional in Eg. (6), the following solution can be 
obtained: 
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Oe a ee ee eee 
i”) , 2, © {24 
t2 ¢3 5¢ 
X, th=t-- +—+—.., (22) 
at=t-5 +E tag 
to t4 
X,t)=—-—+ 
3 t) eee 


4 | Conclusion 


In this study, the differential transform method is successfully expanded for the solution of chemical 
kinetics problems. Since the DTM gives rapidly converging series solutions, the differential transform 
method is more effective than other methods. Finally, all results prove the validity and efficiency of these 
methods in solving nonlinear differential equations. 
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